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Abstract 

The free Schrodinger equation has constant velocity wavepacket solutions ■i/'v of the form ■0 — /(r — 
vi)e-*™'='*/2. These solutions are eigenvectors of a momentum operator p which is symmetric in a pos- 
itive definite scalar product space. We discuss whether these t/^v can act as basis states rather than the usual 
plane wave solutions. 



1 Introduction 

We establish a class of solutions to the Schrodinger equation that are localised, their amplitude inversely 
proportional to the distance from the wavepacket centre; they are non-dispersive, and are eigenfunctions of 
a particular momentum operator found below. We show that the Schrodinger equation 

{idt -I- VV2m) V = (1) 

(we put h = I throughout ) admits the following bounded solution representing a wavepacket travelling with 
velocity v in the z direction 

sin (m{c^ + v^y/^ (x^ + + (z - vt^y/A ^ 

j V /_ i m V [z — vt) — imc t/2 

or more generally with velocity v : 

_ Sin(mt;°|r - V^|) ^„„v.(r-vt) „-^mcH/2 

^v- |^_^^| e e {1} 

with 

i;" = (c2 + |vni/2^ (c2 + „2)i/2^ r = r-vt, f = |r| , 

and jo is the spherical Bessel function of order zero. The c is an arbitrary parameter having dimension 
velocity: it can take any finite value or be zero, as in an earlier version of this paper. From now on we will 
consider c to the speed of light then 

mcv^ — mc (c^ + v^y^"^ ~ m(? -f for w ^ c , 

which is the non-relativistic energy after subtracting out the 'rest energy' rac^ . In the next section we find 
an operator fp such that p^V'v = ip^ . 



To verify that (2) is a solution of (1), note the identities 

"sin(mti''|r — wt 



2 

V +-^)Vv = (3) 



L |r-vt| 
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The first identity (3) follows from ■^v = /(r — vf) e . Expanding out (4) and substituting in (3) we 

obtain 

(s/'^ - 2imv- V - m^v^^il)^ = -m^(c^+t;^)Vv 

2 

{v"^ + 2im{dt + ^"Y") ) V'v = - m^c^ Vv 
(V^ + 2imd^ Vv = 

which is (1). 

That V'v = /(r ~ vt) e~*™'=^*/2 is a clear statement of its non-dispersive property, the e~""'=^*/-^ being 
merely a time phase factor over all space. The wavepacket centre is at r = . The only free parameter in ■0v 
is V itself (once we have fixed the value of c ) , which makes us curious as to whether is an eigenfunction. 
We will go on to consider whether we can regard these Vv as the basis states of the Schrodinger equation, 
rather than the usual plane wave solutions which are spread out over all space. 

I have not found any reference to the V-'v solutions in the literature: it is sometimes stated that such 
constant velocity non-dispersive wavepackets do not exist (see for example pl67 of ref [1]). Wavepacket 
solutions of the Schrodinger equation are considered in Besieris et al [2] (see (2.7) therein), and Barut [3]. 
A constant acceleration wavepacket [4], sometimes called the Airy packet, is well known. 



2 The velocity (momentum) operators 

The momentum of a particle is simply the product of its mass and velocity p = mv as usual for the 
non-relativistic case. As u*^ = \/c^ + we will label 



The has relativistic character, even though the absolute velocity has no upper limit. In terms of p the 
wavef unction (2) is 

_ sin(p°|r - pf/ml) ^^p.^i-.p^/^^ „-»mc^t/2 
|r — pf/TO| 

and at time t = 

V-v ^ - sin(/r)e'P-'' = — - L^p" '■+^p-- - e-'P°'-+'P--l . (5) 
r 2 r L J 

We first establish the following identities, which can be verified by direct calculation: 

i 'r(-m2c2 + V2)(e±'P°VP--) = ± i (a,r)/ (e^'^" e'P-") (6) 
^^giipVgip.r^ = (ip ± i/f) (e±^P°''e'P'") 
(a^r) V(e±'f°''e*P-'') = (zp ± zp"f) (a^r) (e±'P° e'P"') (7) 

now we multiply (6) by f = r/r and subtract (7) obtaining 

[-(a^r)V - ir (m^c^- V^)] (gi^p" p"") = - « (a,.r) p (e±*P" V;'P"-) 

[-(a^r)V - \v(m^c^- V2)]{sin(/r)e*P-'-} = - i (a^r) p {sin(/ r) e^P'''} 

(- \ rr?(?v + a) r V'v = - « {prT) rpip^ . (8) 

The operators a, and aP with (a°)^ = a^, are 

(a°,a) = (-irV^, -(a.r)V + irV^). (9) 

The operator a is related to the Runge-Lenz operator used for solving the Schrodinger equation with a 
Coulomb potential, its components commute with each other. Further properties of (o°, a) are listed in [5]. 
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We now introduce the dilation operator S 
T,= -i (drr) 

and an inverse dilation operator : 

j:-'f{r,9,ct>) = -[ f{t,e,4')dt, E-V(r)=i/ /(Ar)dA. (10) 

1" Jo Jo 

The operators S, commute with any operator of homogeneity degree zero, i.e. any operator such as 
(rV) which is invariant under a dilation of r . Now we can write (8) as 

(— ^m^c^r + a) rV'v = Srp^'v 

and multiplying from the left by ^ we obtain 

pV'v = 7^-1 (- im^c^r + a)rV'v = pV'v, 

p = i5]-i(-imVr + a)r (11) 
which is the momentum operator. 

The operator 

We look for an operator ^ such that p° = P° V'v ? a more difficult task than finding p . As the notation 
implies, we will show in the next section that (p°, p) is a 4- vector, finding the boost operator generating a 
change in momentum. Prom (6) we find that 

(|mVr + aO){sin(A)e'P''} = ^rirn^c^- V^) {sin(A) e'P"'} = - (arr)p°{cos(p°r) e'P-""} 

(12) 

which cannot immediately be resolved into an eigenvalue equation due to the cosine function on the RHS 
instead of a sine function. In Ref [6] we constructed a Hilbert transform operator H- which converts the 
{cos(p° r) e'P""} into {sin(p° r) e'P""} . This H- is defined 

W±/(r) =^[{ne + Ho) ± {He -no)V] /(r) (13) 



where He, Ti-o are the Hilbert transforms of even, odd functions: 

Note that 'H± involves an integration over the entire axis r = A f , (— oo < A < oo) through the origin to the 
field point. We list some properties of these operators (for details see [6]) to be used later, which may be 
verified from the definitions (13) and (14): 

TieTi-o = 'Ho'He = ~ 1 =^ Ti-Ti+ = Ti+Ti^ = — 1 , (15) 
drH± = H^dr, n±r = rH^ n±r = rH^ (16) 

VH±=H±V, n±r = rn± 

il-He/or)^ = -iHo/er {^n±r)^ = -^n^r. (17) 

(15) and (17) imply that {i - r) is a unitary operator. 
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Returning to the operator we first veriiy that W+{cos(p'^ r) e'P ''} = — {sin(p'^ r) e'P ''} , we will use 
Ho sin(Aa;) = cos(Aa;) , He cos(Aa;) = — sin(Aa;) for A > , also noting that {p^r + p • r) > . 



W+{cos(A)e'P--} = - [{He + Ho) + {He-Ho)V] {cos(pOr)e'P--} 

= i [He{l+V) + Ho{l-V)] {e'(f°'' + P--) + e-*(f°''-P--)} 

= i He{cos{p^r + p • r) + cos(p°r — p • r)} + i Wo{sin(p°r + p • r) — sin(p°r — p • r)} 

= ^ — {sin(p"r + p • r) + sin(p'^r — p • r)} + i {cos(p°r + p • r) — cos(p*'r — p • r)} 



^i(p°r + p-r) _ g-i(p°r - p-r) 

sin(/r) e'P'') 



and conversely 

H_{sin(p°r)e'P "'} = {cos(p°r) e'P-"} . 

We substitute this last result into (12) obtaining 

(imVr + aO){sin(A)e'P-''} = - (a^r)p°W_{sin(p°r) e'P""} 

now multiply (20) from the left by — « ^ 7Y+S^^ to obtain 



(18) 
(19) 

(20) 



so that 



pO= -aw+S-i(imVr + aO)r 



1 i 

2 r 



W+S-V(m^c^ - V^); 



(21) 



and we have now found the complete set of velocity operators (p*^, p) . Note that the operator on the LHS of 
(20); which in \r (m^c^ — V^) r , is by inspection a positive operator; and we will later show that (— drvH-r) 
on the RHS of (20) is also positive. The operators (11) and (21) are the energy-momentum 4- vector operator 
p^ : 

= (/, p) = (- iiW+S-i(im2c2r + aO)r, i S"! (- i m^c^r + a) r ) (22) 
p'^ = (a/ m^c^ + p) 

As we recall from (2), the eigenfunctions of contain the same plane wave eigenfunction e*P '' of the usual 
momentum operator — iV , except that our ip^ also has the envelope factor .7o(j'"^) • We next find a set of 
boost operators K and rotation operators J which together with p^ satisfy the Lorentz group commutation 
relations. 
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3 The Lorentz operators 



We establish the Lorentz group operators which must obey the commutation relations 

[J", p^] = i e^^^p" , [K", fP] = ip", [K\ = i (J'^y , (23) 

[J°, Jb] = i e'^b'^J" , [J'^, K^] = i ("'"'K" , [K", K^]= -i e^^^J" (24) 

with the antisymmetric tensor e^^^ = 1 . The boost, rotation operators K, J such that p^ is a 4-vector are 

K=-VH-r, J=-i(rxV), (26) 

the rotation operator is as usual, but the boost operator necessarily contains the 7i_ integral operator. We 
first verify that the components of K satisfy the last formula of (24), using the identities (16) 

[K^^ K2] = [- ViH_r, - VaH-r] = ViH-rVaW-r - VaW-rViW.r 

= ViH-rH-Var - VzH-rH-Vir = ViH-H+rV2r - WiU-H+rWir 

= - VirV2r + V2rVir = - (r x V)3 = - iJs 

as stated. 

To verify (23) we first separate out the terms containing in the (p°, p) operator of (22) as follows: 



(/,p) = iS-VmV(iW_r,r) + iS-V ( - iW+a° r , iar). 



(27) 



The Lorentz operators (26) commute with ^Sr , so we now have to verify that both 

(iW_r, r) and ( - iW_a° r , i ar ) 

are 4-vectors. In the first case we have 

[K, iH-r] = [-VH-r, iH-r] =-iVn-r'H-r + iU-rVH-r 

= -iVrH+H-r + irVH+H-r = iVrr - irV r = ir (28) 

and [K", = -V'"H-rr'' +r''V''H-r = -V''r'"H-r +r''V'"H-r = iS^^iiU-r) 

as stated. This result indicates that r is the space part of a 4-vector: the space-time implications of this we 
will discuss elsewhere. Also 

[K, -^n+a°r] = [- VW_r, -^H+a^r] = iVH-H+a'^r - i^H+a^rVH-r 

= - i V a° r - i ^H+H-oP rV r = -iV aP r + i ^ rV r (29) 

= i[- irV , o°] r = i ^ ar , 
the commutator on the last line is a known result (see for example [8]). To derive (28,29) we have used (16): 

drn± = n^dr, n±r = ru^ n±r = rn^, vn± = n±v, n±r = rn±. 

We next investigate the effect of the boost operator K of (26) on tp^ ■ Under an infinitessimal boost e 
in the z direction the wavefunction at time i = is transformed as [V'v]j_q ^ (1 + i^K^) ['0v]j_q- We 
use (19) to obtain 

"sin(p°r) 



. N r , n sin(B°r) . , „ „ , 

(1 + leKz) [V'v]t=o = — exp[ip-r] - leVsH-r 



ex.p[i p • rj 



sm{p°r) 



exp[i p • r] — i e V3 cos(p°r) exp[i p • r] 



/, • n X sin(p"r) 
(1 + lep^'z) — — - 



ep3 cos(p°r) 



exp[i mv • rJ 



We can show that this is equivalent to substituting 



P3 



P3 + ep , Pi^Pi, P2^P2, P -^P +ep3 



(30) 
(31) 
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into Vv , as then 

r , 1 sin((p° + eps) r) q 



sinfo^r) cos(p'^r) 



r r 



exp[i ep°2:] exp[i p • r] 



• n ^ Sin(p°r) / n ^l r- 1 

(1 + lep z) h ep3 cosm r) expkp • rl 

r J 



(32) 



which is (30). This confirms that K generates the transformation (31). 



4 The positive definite scalar product space {(f) | tp)g 

We estabhsh in the Appendix that 

{—drV Ti- r ) = {—i S Ti- r) is a positive operator. (33) 

Recalhng (17): 

(-H±r)t = --n^r, 
r r 

and also 

r 

then 

(-iSW_r)t = (-iSr(iW_r))^ = (-i(iW+r)Sr) = (-iS(iW+r)r) = (-iSW_r), 

so that {—i S W_ r ) is self-adjoint as well as positive. So we can construct a positive definite scalar product 
space, a Hilbert space: 

((/)| V)s = -i^H-ri)). (34) 

We will verify that the K, operators of the last section are self-adjoint with respect to | ip) ^ . (The 
rotation operators commute with the (— iEH_r) operator). We show the working involved to check the 
adjoint of K = — V 7Y_ r : 

= {(t)\in-Y.rVn-r^) = {- ^^H+r^ (l)\iY.rV H-r^jj) 

= { - ^n+r^ (l)\irVY.'H-r-4)) = ( V ^H+r^ ^ | i S H- r V) 

= {V'H-r<t)\iT.'H-rxl}) = (K^|^/;)g. 

In going to the last line of working we have used T-L± r = rHzf from (16). In a similar fashion the self- 
adjointness of the can be verified with the adjoint property 

(a^r)t = (a^r) . (35) 
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5 The inner product spaces (^(j)\ip) g and | V') • 



We have established in section 2 that the V'v are eigenvectors of which is symmetric with respect to 
the non-negative scalar product space (</> | . So two wavepackets Vv, V'v' momentarily coinciding are 
necessarily orthogonal: 



(V'v I V'v) c 







for V ^ v' 



In the Appendix we check the above result and calculate (V'-v | V'v)^ and (Vv | V'v)g • We find that 

(^v|^v')5= (2 7r)3^5(p-p') = 47rV<5(p-p')- (33) 



The physical meaning of the scalar product space (V" | V')5 is obscure. More interesting is the usual scalar 
product of two wavcfunctions (V | V') ) the {(fftjj) within the integral being a conserved quantity for any 0, if^ 
satisfying (1). Consider two wavepackets ipyi with velocities v, v' in the z direction which momentarily 
coincide at time t = , so that ipy = sin(p°r) e'^^/r and V'l;' = sin(p'° r) e'^ ^/r , with — oo < {p,p') < oo . 
Then 

<V'^'|V'^)= J {sm{p'\)e-'P'yr){sin{p°r)e'Pyr)d'^r 

= 27r y sin(p'°r) sin(p°r) (f-f ^ d{- cos 6) dr 

sm{{p - p')r) -\ 
{p -p')r J 

" sin([p -p' + (pO - p'0)]r) + sin([p - p' - {p° - p"')]r) 



cos((p° -p'°)r) - cos((/ + p'°)r) 

- sin([p° +p'° +p- p']r) + sm{[p° + p'° - p + p']r) 
sgn(p - p') + sgn(p - p') - 1 + 1 



- dr 
r 



ip-p') 2 

^2 



" \{P-P')\ ■ ^^^^ 

So that the scalar product of two wavepackets which momentarily coincide have an inner product (Vt^' | V'v) 
which is a finite number (when v ^ v'), and furthermore this inner product is constant for all time, however 
far apart the wavepackets are separated (either before or after the wavepackets have coincided). This linkage 
between two coinciding wavepackets or particles is a well known quantum feature. For two wavepackets 
coinciding with momenta p, p' not necessarily collinear, the above result (34) is 



(V'v I V'v) 



lp-p'l 



(35) 
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6 Towards a fully relativistic theory 



The evolution equation for the ip^ is the non-relativistic Schrodinger equation, and the momentum p = mv 
which imphes an unbounded velocity. On the other hand the energy cp° — ca/ m^c^ + as in the relativistic 
case. It is straightforward to obtain a so that the 4-momentum has the relativistic value 

p) *v = {mjc, m7v) , 7 = (1 - vVc^y^^^ , 

but more difficult to find the evolution equation for this 

Consider 

= — ^ ' exp[zm7v • (r - vt)] e ^'^^ (36) 

then from the results of section 2 by simple substitution 

(P° , P) [*v] = (™7C, m7v ) [*v] (37) 

and now the (while still having velocity v ) has the correct relativistic 4-momentum when acted on by 
the operator. We look for the evolution equation for • Proceeding as in section 1, we note the identities 

2 

(a^+vV = (38) 

sin(m7|r — wt\) 



(V - im7v)«'v = e*'"^-('--^*)e-'™^'*/2'y V ^ _ 

(V - im^wf^^ = e'™v.(r-vt)g-imc=t/2^ ^2 [ sm{m^\r - vt\) i ^ 2^2^^ (39) 

r — 



Expanding out (39) and substituting in (38) we obtain 

(s/^ - 2im7v- V - m^7^t;^^ ^'v = -rn^7^c^*v 

9 

(1 Tnc \ 
+ 2im-f{dt + -^^)) *v = -rn^c^^v 

(V^ + 2im-fdt^ *v = (V^ + 2i^at) *v = (40) 

which apart the gamma factor in the partial time derivative is just the Schrodinger equation (1). But due 
to the 7 factor which depends on the wavepacket velocity, equation (40) is not a true evolution equation 
applicable to any ^'v Wc can try inserting the p° operator instead of p° in (40), but the jP operator only has 
the eigenvalue property when the origin is at the wavepacket centre, which is not desirable for a evolution 
equation applicable wherever the origin. 
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7 Outlook 



The attraction of the wavepackct solutions ip^ is their locahsed nature which accords with our everyday 
experience. We have made progress in showing that the Vv can be used as alternative basis states to the 
plane waves, albeit with some complications: for example the momentum operators (22) only have the 
eigenfunction property when the origin is at the wavepacket centre. 

The usual momentum operator tt = — i V does not commute with the momentum operators p: the 
wavepackets are not invariant under translation, unlike a plane wave - which merely undergoes a phase 
change on translation. (Of course the Vv are invariant under a combined space and time translation.) The 

■i/)v arc not cigcnfunctions of the Hamiltonian operator — V^/2m : it is impossible to have stationary states 
of the form %jj = e^*^*/(r) for wavepackets in motion, as the wavepacket envelope must itself be a function 
of t . An unexpected feature of the wavepacket momentum operator is that p is part of a 4- vector with 
m^c^ + (p)^ = {p^Y ) as discussed in section 3. From (30) we see how the boost operator generates a 
change in momentum. However the velocity has no upper limit, the operator p/m is unbounded. In a 
fully relativistic theory we would expect the wavepacket velocity to be less than c ; so far, as shown in the 
previous section, we have had limited success in the relativistic wavepacket extension. We have shown that 
the position r is part of a 4- vector with r° = z Ti-r , which has far-reaching implications on the nature of 
space-time. We will discuss this further elsewhere. 

Operators very similar to p have been discussed before in the context of lightcone quantum mechanics, 

an approach initiated by Dirac [10] in 1949 and subsequently worked on by a number of authors [7,8,9]. 

Their approach was to start with the Poisson bracket formulation of classical mechanics, find the lightcone 

generators, and then 'quantise' these generators to appropriate operators. The second of these authors [8] 

found a canonical transformation from the usual conjugate pair (r, tt) (x, p) , where 

, , mV-FTT^ , , 2(r-7r) 

p(r , tt) = TT — — r , x(r , tt) = „ , , — ^ tt . 

2 (r • tt) m-^c^ + tt^ 

Our p operator appears to be the quantised version of the generator above with tt ^ — i V , which is more 
easily seen if we write p in the form 

p=-ziVr-iiE-ir (w^c^ - V^) r , 
the S corresponding to the classical (r • tt) . Consider a positive energy plane wave g-^cp^t+ip r ^ q^-^ ^j^g ^ast 
lightcone t = —r/c this becomes e*P°''+'P''' , which referring to (5) is one half of the ip^ wavefunction, the 
other half containing e~*p'''"+*P''' can be thought of as a future lightcone wave. Prom this point of view the 
V'v is a balance of past and future lightcone waves. The wavepacket approach hints at a more dynamical 
picture of a particle - existing as a confluence of past and future lightcone waves. 

Acknowledgement I am grateful to Prof G Kaiser [11] for pointing out that puting the light velocity 
c = 1 as in an earlier version of this paper leads to incorrect dimension. 
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Appendix 

I. The positive operator {—drrTi-r) 

That {—drvH- r ) is a positive operator is a key result in this paper, leading to the Hilbert space of section 
4. We reproduce the derivation from ref [6] below. 

Consider the unitary operator U 

Uf{r) =^[i^c -iJ's) + {J'c + i:F,)V]r f{v) {Al) 
where V is the parity operator and J^c, are the Fourier cosine, sine transforms defined by 

P/(r,^,(/.) = /(r,7r-^,(/. + 7r) 



f{r, e, (t>) = 9c{r, 0, </.) = W - / f{t, 9, </.) cos(ri) dt , 



Ts f{r, e, 4>) = 9s{r, 6,4>)= ^ f{t, 6, 4>) sin(rt) dt . 
The operator U has inverse where 

U-^ = U* =^ [{:Fe + iJ's) + {To - iTs)V\r. (A3) 
It is straightforward to check that U*U = 1 , with the aid of 

TcTo=TsTs=VV =1. (A4) 
We emphasise that lA^ U* are in a sense one dimensional transform operators: the integration is along the 
axis Ar , with — oo > A > oo . 

We next calculate the operator 

f=WrU (A5) 

which is clearly a non-negative operator. To simplify (A5), it is known that the cosine transform of the sine 
transform is the Hilbert transform of odd functions, and the sine transform of the cosine transform is minus 
the Hilbert transform of even functions, i.e. 

TsTc = — We , TcTs = 'Ho 

where He, Ho are defined by 



f( a,.^- r IMiAm f( ^ 2 /(Ar) 

n^^- 2 r tf{t,e,ct>) ^ ^, , 2 pA/^Ar) 



For brevity we adopt the notation 

T± = Tc ± iTs 

then from (A4,A6) there follows the further identities 

T+T+= -i{He-Ho), T-T-^ i{He-Ho), 

T+T- = 2 - iiTie + Ho), T-T+ = 2 + i{ne + no) 



Also we note that 

T± r = Tidr T± , r T± = ± i T± dr ■ (A9) 

Returning to f = U* rU we have 

=U*rU =^ {T+ + T-V )r {T- + T+V ) r 

= i-^dr + T-V) {T- + T+V) r 

4r 

= i^dr [{-T+T- + T-T+) + {T-T- - T+T+)V\ r 
4r 

= -^Qr [(He + Ho) + {He-Ho)V] T = --OrH+r {AW) 
2r r 
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recalling W+ = i [(We + T^o) + (We — Wo)^] from (13) which is a Hilbert transform along the axis 
Af , (—00 < A < 00) . 

The operator components dr , W+ of Z° do not commute: If in the working out of (AlO) we take r to 
the right instead of to the left, we arrive at 



f = U*rU = - ^ [{Ue+Uo) - {ne-no)V] drr 

= --H-drT. 

r 



As 



1, 



1 , then 



which can be verified their definitions (13) together with TieTi-o = Ti-oTi-e - 

= -drU+H-drr = --dlr 
r r 

as expected. Hence /° is a positive square root of —^d^r. 
Finally from (AlO), (All) 

rfr = {-drH+r'^) = {-drrH-r) 
and as we have shown 1° is non-negative, so must be the operator rl^r = {—drvH- r ) . 



(All) 



{A12) 



II. To calculate (V'-v | V'v)^ 

We calculate (^tp^ | ipv')g, the inner product of two wavepackets coinciding momentarily at time t = 0. 

If V ^ v' then there will be a centre-of-momentum frame where v = — v' . (The exponentiated boost 

operator e*^'^ can be applied to the V'v, V'v' to arrive at the centre-of-momentum frame.) Before we calculate 

{t/j-v \'4'v)g we first simplify the inner product (0 | V')^ : 

{(f)\tlj)s = {(i)\{-iY.H-r)il}) = (0| -dr-rH-ri}) 

= \{ldycl>\H-ri,) + \{cl>\ -drrH-ri>) 

= \{{dr+ i)r<^|W_rV) + \{<l>\- r{dr+ i)W_rV) 



= \{dr{r4>)\n-ri,) + \{r<j>\- dr{H-r^)) 



(A13) 



(9, (r </)*)) (H_rV) + (rr)(-a,(W_rV))^ 

Now we substitute if) = ip^, 4> = V'-v into the above, recalling from (19) that H-ripv = {cos(p°r) e*P ''} , 
obtaining 



(^--v I '0v)< 



dr ( sin(p°r) p") (cos(/r) e' p") 

- ( sin(p°r) e' p") {dr cos(/r) e' p") 

-— / (cos(p°r) -|- ■ f) sin(p'^r)) cos(p'^r) 

■ I P 



+ sin(p°r) (sin(p°r) — % (— • f) cos(p°r)) 



'2 



j e'^'P-'Uh = (27r)^^5(p- (-p)) =0 for 0. 



{AU) 



If instead we substitute ip = (f) = ipv into (30) we obtain 

P^ f \ P 
(V'v I V'v)_5 = y y (cos(p°r) - i (^ • f) sin(p°r)) cos(p°r) 

+ sin(p°r) (sin(p°r) — i (^ • r) cos(j3°r)) 



d-'r 



^ y [1 - 2 i ( J • r) sin(pOr) cos(/r) ] d^v = (2 Tr)^ ^ <5(p - p)) , {AW) 



11 



as integrating over the angular variable (p • f ) = cos 6 annihilates the second term in the integral. We 
conclude that 

(^,|V-v')s= (2 7r)3^^(p-p') = 4^V'5(p-p')- (^16) 
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